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Abstract 

Several series of shallow large-size (Z 3> lattice spacing) two-hole bound states 
are found in the two-dimensional t — J model. Their binding energies depend ex- 
ponentially on the inverse value of the hole-spin- wave coupling constant. Their 
parity is connected with the angular wave function in a nonstandard way. The 
possible role of these states in the formation of superconducting pairing is con- 
sidered. 
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I. INTRODUCTION 



The problem of mobile holes in the t-J model is closely related to high-T c superconductivity. 
The t-J model with less than half-filling is defined by the Hamiltonian 

H = H t + Hj = -t £ (dl a d ma + R.c.) + J £ S n S m , (1) 

<nm>cr <nm> 

where d^ is the creation operator of a hole with spin a {a =f, J,) at site n on a two-dimensional 
square lattice. The d\ a operator acts in the Hilbert space where there is no double electron 
occupancy. The spin operator is S ra = ^S na a a ^d n ^. < nm > are neighbour sites on the lattice. 
Below we set J = 1. 

At half-filling (one hole per site) the t-J model is equivalent to the Heisenberg antiferro- 
magnet model which has the long-range antiferromagnetic order in the ground state 0-§J. 
The problem is the behavior of the system under doping by additional holes. 

Single particle properties in the t-J model are by now well established. A single hole is an 
object with complex structure due to virtual admixture of spin excitations. It has been shown 
that a single hole has the ground state with a momentum of k = (±7r/2, ±7r/2). The energy is 
almost degenerate along the line cos k x + cos k y = which is the edge of the magnetic Brillouin 
zone (see e.g. Refs. [§-|T3|). 

The two-hole problem is much more complicated. There is no doubt that for tcl there are 
short-range d- and p-wave two-hole bound states with a binding energy AE w —0.25. These 
bound states could be relevant to the superconductivity. However the value of t corresponding 
to the realistic high-T c superconductors is t ~ 3 (see, e.g. Refs. Hl6-18fl). The early works on 



exact diagonalizations for finite-size clusters [|ToyT9|-[23| indicated the existence of a short-range 



two-hole bound state for large t. However recent investigations ||2~4H27| have demonstrated that 



the short-range two-hole bound state vanishes at t w 2 — 3. Moreover, in our opinion even for 
smaller values of t the short-range Coulomb repulsion in realistic systems should destroy this 



bound state (see the discussion in Ref. f27j). These arguments lead us to the conclusion that 
theory predicts no short-range bound states at physical values of t. This conclusion makes the 
investigation of the long-range dynamics in the t — J model very important. 

The long-range dynamics is determined by the interaction of a hole with spin-wave exci- 
tations. Calculation of the hole-spin-wave coupling constant / is straightforward if one treats 
H t as perturbation and uses the condition zS ^> 1 (z is the number of neighbour sites). The 
value of the coupling constant obtained in this way is f — \zt = 2t (see e.g. Ref. P, [H2| , |r5fl ) . The 



calculation of / for realistic zS and arbitrary t has been done in recent works [25LBD| 



In the present work we prove the existence of shallow large-size two-hole bound states which 
appear due to the exchange of a spin-wave. It is demonstrated that there are several infinite 
series of such states. Their binding energies depend exponentially on the radial quantum 
number: E n oc e _77rn , where 7 is a numerical constant, n = 0, 1, . . .. Very interesting is the 
relation between the parity and the angular wave function of these states. A state with an 
angular wave function without nodes (in this sense it is s-wave) has negative parity. A state 
with two nodes in the angular wave function (in this sence it is p-wave) has positive parity. 
This property may be considered as an additional negative "internal" parity of the pair. The 
nonstandard parity relation is caused by two physical reasons. First, the ground state of the 
system is the antiferromagnetic one. Second, the bottom of single-hole dispersion is at the face 
of the Brillouin zone. 

Our paper has the following structure. In Sec. II we derive an effective Hamiltonian describ- 
ing the dynamics of the system in the long-range region. It is presented in terms of dressed 
holes and spin-waves. The nontrivial properties of the system are caused by the strong interac- 
tion between holes and spin-waves and the contact interaction between the holes. This section 



is based on the results of Refs. p7H29|]. In Sec. Ill we show that the exchange of a spin- wave 



between two holes creates the long-range attracting potential which decreases as 1/r with the 



separation r between the holes. This potential generates several infinitive series of bound states 
of the two holes. In Sec. IV we discuss the parity and the symmetry properties of these states. 
In Sec.V we compare the analytical solution with the results of exact numerical diagonaliza- 
tion of the equation for the bound states. In concluding Sec. VI we discuss scenario of possible 
superconducting pairing based upon the results of the present work. 

II. INTERACTION BETWEEN THE HOLES WITH OPPOSITE SPINS. 
EFFECTIVE HAMILTONIAN FOR THE DRESSED HOLES. 

We treat the J term of the Hamiltonian ([!]) using the linear-spin-wave approximation (see 
Ref. [Bfl for a review). Define the following Fourier transformations 



' net ' nel 



where our notation n ef (n ej,) means that site n is on the up sublattice (down sub lattice). 
Introducing the Bogoliubov canonical transformation, 

«q = « q a q - v q 6_ q , /? q = w q 6 q - v q a_ q , (3) 

we write the Heisenberg Hamiltonian Hj as 

Hj = E + J2^(^ + 13%), (4) 
q 

where E is the antiferromagnetic background energy. The summation over q is restricted 
inside the Brillouin zone of one sublattice where 7 q = |(cosg x + cosq y ) > 0. The spin-wave 
dispersion and the transformation coefficients are given by 



cu q = 2Jl - 72 -> v^lql, at q < 1 



u n = \ — + (5) 



w q = -sgn(y, 



2 



The spin waves a q and /3 q have definite values of the spin projections. Due to Eqs.(|]),@ acUO) 
has S z = — 1, and Pl\0) has S z = +1. Here |0) is the quantum Neel state wave function. 

A single hole on the antiferromagnetic background is a composite object due to strong 
virtual admixture of the spin excitations. Therefore its dispersion is sufficiently complicated. 
Fortunately it may be well approximated by the simple analytical formula which was found in 



Refs. MM. 



e k = ^A 2 /4 + 4t 2 (l + y) - ^A 2 /4 + 4t 2 (l + y) - 4t 2 (x + y) 7 2 + -/3 2 (cosk x -cos k y f (6) 
A « 1.33, x « 0.56, y « 0.14 

The parameters A,x,y are some combinations of the ground state spin correlators The 



fa term gives the dispersion along the edge of the Brillouin zone. It is rather weak and may be 
estimated as fa ~ 0.1 x t at t > A/4. Near the band bottom k = (±7r/2, ±vt/2) the dispersion 
(H) can be presented in a usual quadratic form 

e k « -ftfc? + ^ 2 2 , /3 2 «/3i, (7) 

where &2 is projection of the momentum along the face of magnetic Brillouin zone, k\ is or- 
thogonal projection of the momentum, see Fig.l. Due to Eq.(|5]) 

4^t 2 « 2.lt 2 for t < A/4 = 0.33, 

(8) 

-#t«0.65t for t> A/4. 
A comparison with numerical calculations (see e.g. Refs p.pd ^l^IT5|j29|| ) shows that simple 



formulas ©,@ underestimate Pi by 10-20%. This produces no serious effect on the results of 
the present work, but certainly this correction can to be included. 

The wave function of a single hole may be written in the form ip^ = /ij^jO). is the 
creation operator of a dressed hole 

hU = Mj2hl^, (9) 

' nej. 



hi 



' net 



It was demonstrated in Refs. |15| , [29|] that may be well approximated by the following simple 



expressions: 



V 



3/2 + 2S t 
S t 



1/2 



[S t (3/2 + 2S t )]V 



2 ' 



(10) 



5 t = v/9/16 + 4t 2 , 



where 5 is a unit vector corresponding to one step in the lattice. This form of trial wave function 
leads to the dispersion (^]). 

The effective interaction of a composite hole with a long wave-length spin-wave is of the 
form (see, e.g. Refs. ||00[||[29l) 



9(K q) ( + 4+qT^kl/^q + H.C. I , (11) 

k,q V ' 

g(K q) = 2/^/^( 7k U q + 7k+q^q)- 

In the perturbation theory limit it <C A/4) and at zS 3> 1 [z is the number of neighbour 
sites) the coupling constant is / ~ f = 2t (see e.g. Ref. P |L2T|L3"1| ). Account of the first 1/zS 
correction gives / « f% — |t, and summation of all 1/zS series gives / = 3.4t for f <C A/4 



. For an arbitrary t the coupling constant was calculated in Ref. [29|. The plot of / as a 
function of t is presented in Fig. 2. For small t: / « 3.4t. For large t the coupling constant is 
t-independent: / ~ 2. 

The interaction between the two holes can be caused by the exchange of single (or several) 
soft spin-wave. Alongside with that there is a contact hole-hole interaction. One can say that 
it is due to the exchange of several hard spin-wave excitations. Now we are going to derive 



the Hamiltonian of the contact hole-hole interaction. Interaction between two holes at large 



momentum transfer (q ~ 1) has been obtained in Ref. [27| using a variational approach and 
1/zS expansion. 



r(k 3 T,k 4 |;k! T,k 2 



c 

^Tid-ks + 5(7i7s + 7274) + ^-(7k 1+ k 3 + 7k 2 +k 4 ) 



; 12,34, 



'121 



A = 16tz//i 3 (l - 7^ 2 ) - ~ - 2fi 2 - 18.5/i 4 + 84/i 6 + lOatz/ 3 // 3 , 

B = 8tz//i(l - 9/i 2 + 32/i 4 ) - 4/i 2 (l - 4/i 2 )(l - 12/i 2 ), 

2 ~ 
C = -atv\i . 

Here v and /i are the parameters of the wave function (|i~0"[). We denote for simplicity 7^ = ji- 
The functions A and C in formula fli~2]) have been derived in the first order in a, where a 
is the coefficient in front of the transverse contribution to the Heisenberg energy: S n S m — > 
Sn^m + f {^n^m + ^n^m)- The physical value is a = 1. The higher order contributions to A 
and C are essential. In order to fit them one has to set a in Eq.flT2|) somewhere in between 
a = 0.5 and a = 1. The expression for B in formula (|1^) corresponds to the Ising background 
(a = 0). The difference is that A(a = 0) and C(a = 0) are small and therefore the corrections 
are crucial. The coefficient B is not small and the a correction is not so important. 

Expression ([12]) includes the contributions from the single spin-wave exchange as well as the 
many spin-wave exchange. It is important to separate them because, as it will be clear from 
further consideration, they play quite a different role in the physics of large-size bound states. 

Let us demonstrate now that the B-term in (|T2|) is due to the single spin-wave exchange, 
see Fig. 3. For the Ising background there is no hole dispersion, the frequency of the spin-wave 
is equal u; q = 2, and the Bogoliubov's parameters in Eq.(||) are w q = 1, t> q = 0. Therefore due 
to Eq.(Jll"|) the single spin- wave exchange contribution (Fig.3) is equal 

8/ 2 1 4/ 2 

- ~TT~/ 7(7i73 + 7274) = -rr (71 73 + 7274)- (13) 

N (-o> q ) N 

We see that the single spin-wave exchange (0) possesses the same kinematic structure as the 
5-term in (|T2"|). Let us verify that their absolute values are also identical. The interaction ([12]) 



has been derived in first order in 1/zS expansion. In this approximation for t « A/4 the hole 
spin- wave coupling constant is / = |t (see Ref. [p8|). Substituting this value into expression 
(13) one can easily check that it perfectly agrees with B = ^-t 2 obtained from Eq. (fl~2|) . At 
large t the agreement between (|T3"|) and 5-term in Eq. (|T2|) cannot be perfect because the latter 
comes from an approximate variational solution. Nevertheless numerically the agreement is 
good enough. For example for t = 3 the coupling constant is / ~ 1.8 (see Fig. 2) and the value 
of B obtained from formula (0) is B « 1.6. From expression (|T2| ) one gets B w 2.2. This 
consideration demonstrates that the 5-term in the short-range interaction ( |T2"| ) accounts for the 
single spin-wave exchange. Therefore A and C terms describe the contact hole-hole interaction. 
Due to formula fO) the Hamiltonian of contact interaction is of the form 



Hi 



I'!' ~ iu E 

1,2,3,4 



c 

^7k!-k 3 + -^ki+k 3 + 7k 2 +kJ 



(14) 



We denote for simplicity = hi. 

Let us formulate the result of the present section. Starting from the Hamiltonian of the t-J 
model (0) which is expressed in terms of operators d na we show that the dynamics of holes on 
the antiferromagnetic background is described by the effective Hamiltonian 



H eff = 



ek^k^k + X^q( a i a q + /^qAj) + H h,sw + H, 



q 



Mi 



(15) 



which is expressed in terms of composite operators and spin-waves a q , /3 q . The interactions 
Hh tSW an d H hh are given by Eqs. (|Tl~l) and (0). 



III. TWO-HOLE BOUND STATES 



Consider a system of two holes with total momentum P = and the projection of spin 
S z — 0. Let be the wave function in the momentum representation describing the state of 
the system with the relative momentum 2k: 



|P = 0) = ^=E^4 T ^kjO). (16) 
viv k 

The function g k satisfies the equation of the type of Bethe-Salpeter one: 

(E - 2e k )g k = £ T(E; k, -k; p, -p)g p . (17) 
p 

Summation is carried out over the magnetic Brillouin zone (7 P > 0) . The kernel of this equation 
can be presented in the form 

r Tsui ^contact- 

The term T sw in (|18| ) is the contribution from the single spin- wave exchange (see Fig. 3) for 
which we find with the help of Eq. (pT|) 

F 16/ 2 (7-k^q + 7p^q) (7-p^q + 7k^q) q\ 

iV E - e p - e k - cj q 

where q = p + k. The term r con f act in Eq.(^) desribes the many spin- wave exchange. From 
Eq.flHD we derive 



^contact jy 



^7k-p + C7k+p 



(20) 



Let us apply Eq.(|l7|) to the problem of the large-size bound state. In this case one can 
restrict the consideration to the two-hole dynamics near the band bottom in one pocket of the 
Brillouin zone, for example near the point ko = (n/2,n/2) (see Fig.l). We have to allow k\ to 
be negative as well as positive. The hole state with negative k\ is outside the Brillouin zone, 
but it is equivalent to one with k' = k — g, which is inside the Brillouin zone. Here g = (n, n) 
is the vector of the inverse magnetic lattice. We will show that the single spin-wave exchange 
r sw ( |19D plays the crucial role in the formation of these states. Therefore we consider it first. 
Later we will take into account the term Y contact (p0|). Note also that the binding energy as well 
as the kinetic energy of the pair will be shown to be small compared to the spin-wave energy. 



Therefore we can neglect E — e p — ek compared to u; q in the denominator in Eq . (|T9|) . As a 
result we obtain 



N q 2 N (k + p) 



1 ~ 1 sw ~ Ar „ 2 — Ar /l. , „\2 ' 



We have shifted the zero of momentum to the center of the poket. Equation fll7|) is reduced 
with the help of fl2lD to the simple form 

(E-Ml-M^Sfj^J^^. (22, 

The last equation is invariant in respect to the transformation k — ► — k, p — ► —p. That is why 
we can classify the wave function by its behaviour under the transformation: 

9k 9-k = Tl9k, (23) 

where we introduce the quantum number 1Z = ±1. Let us stress that 1Z is not the parity of the 
state because k is the deviation from k which is not the center of symmetry of the Brillouin 
zone. The relation between the quantum number TZ and the parity will be discussed below. 
From (H, H) we get 



(E - fck'i - f3 2 k 2 ) 9k = SfU j ^.p^ p- ( 24 ) 



Using the Fourier transformation 

(25) 

we can rewrite it in the coordinate representation as usual Schrodinger equation: 

d 2 „ d 2 



' (3l dx 2 ~ (32 M + u{r)]i) = E ^ m) 



where the potential energy is 



U(r)=n(2f 2 /n)^-^. (27) 



It is convenient to perform the rescaling Xi = y/Pif], £2 — vT^C an d to introduce the polar 
coordinates 77 = pcostp, ( = psin<p. In new variables equation (f26|) has the form 



pdp^dp p 2 dip 2 p 2 ) 



(28) 



where 



. , sin 2 oj — a cos 2 0? „ x 

V(<p) = UFa . J 29 

(sin ip + a cos z if) 2 

F = 2/ 2 /(vr/? 1 ), a = /3 1 //h. 

The plot of the interaction constant F as a function of t is presented in Fig. 2. Typical mass 
ratio in which we are interested is a ~ 3 — 10. The variables in equation (|28|) are separated, 
the solution is of the form 

1> = R{p)*{<p). (30) 



From (28), (|30|) we find the following equations: 

d 2 



dip 2 



+ V(<p))9{<p) = \*{<p), (31) 



1 cf d A 



p- + -)i?(p) = J E;i?(p). (32) 



p dp dp p 

The boundary condition for $(<p) follows from the symmetry relation (^3 



$(<p + 7r) =K$((p). (33) 

The boundary condition for R(p) we will discuss later. 

The crucial thing is the sign of the parameter A, which is the eigenvalue of the angular 



equation (31). The latter has the form of the usual Schrodinger equation in which the function 
V(ip) plays the role of a potential and the eigenvalue A plays the role of an energy. Note that 
V(ip) ( p9"D depends on the symmetry of the wave function (^) through the symmetry parameter 



TZ. If one neglects V(ip) then the solution of Eq.(^TJ) is trivial, A = m? > 0, m = 0, 1, . . .. The 
point is that V(<p) provides the regions of attraction, where V(<p) < 0. For TZ — +1 these 
regions are in the vicinity of points (p = and tp — 7T, e.g. in the direction perpendicular to the 
edge of the Brillouin zone. For TZ — — 1 they are near points tp = tt/2 and <p = — vr/2, e.g. in 
the direction along the edge of the Brillouin zone, see Fig.l. Hence one can expect that there 
are solutions of Eq.(|3T|) with negative eigenvalues A. They provide the attracting potential in 
the radial equation fl3"2| ) and, as a result, the bound states for the two holes. Let us show these 
solutions do exist and examine their properties. 

Consider first the case when the potential V(<p) (^) is small, F < 1. Then the solutions 
of Eqs.([3Tl), (|3lD in zero approximations are &$(tp) = exp(im(p), \$ = m 2 , m = 0, ±1, .... 

For the first level we have m = 0, ^q ^^) = 1, Aq '' = 0. Using Eq . (p9|) one can easily verify 
that first order correction to the eigenvalue vanishes: Aq 1 = J 2?r V{(p)dip = 0. Therefore only 

(2) 

the second order correction, which is definitely negative, is essential: Ao ~ Aq < 0. This proves 
the existence of negative eigenvalue of Eq. (|3T| ) for small F. The considered solution belongs to 
TZ = 1 symmetry. 

Consider now the next level. It is degenerate: m = ±1,$^(<^) = exp(±i<^),A^ = 1. The 
potential V(<p) splits it. The splitting appears in the first order of perturbation theory because 
the level is degenerate. As a result one component of the doublet, Ai, goes down with increasing 
F as: 



Therefore for sufficiently large F one has to expect that this eigenvalue is to become negative. 
The corresponding eigenfunction possesses TZ = — 1 symmetry. 

For large F both eigenvalues Ao and Ai decrease, but Ai does it faster. Therefore for some 
F there is an intersection of Aq with Ai. The reason for this fact is the following. If TZ — — 1 




(34) 



then the potential V(ip) (p9|) has minima near the points tp = ±| which become deeper with 
increase of F. As a result the ground-state function $7£ = _i(<£>) is strongly localized in the 
vicinity of these points. Therefore we can approximate the potential V((p) by the potential of 
a harmonic oscillator near the points ip = This gives the following asymptotic expression 
for the eigenvalue Ai which is valid if Fa 3> 1: 

15a 2 - 15a + 2 



Ax « -Fa + ^Fa(3a - 1) - _ . (35) 

This proves that Ai becomes negative for sufficiently large F. 

The corresponding expression for A is more complicated, because the potential V((p) for 
1Z = 1 can possess either two or four minima for different values of a. Therefore we present 
only the lower estimation for A : 

{ j p 0+ 1 ) 2 for a > 3 
(36) 
—F for a < 3. 

Comparing inequality ( |3~6l) with expression (|35"D we really see that Ai goes down with increase 
of F more rapidly than Ao- 

There is another parameter, the mass ratio a, which governs the value of A. It is clear that 
the higher a the more our two-dimensional system looks like a one-dimensional one. The influ- 
ence of the potential in one-dimension is stronger, the binding energies are deeper. Therefore 
both Ao and Ai decrease with the increase of a. Due to formulas fl3"5]), (35) Ai decreases faster 
than A . 

Fig. 4 presents the results of numerical solution of the eigenvalue problem (|3~ll) , (|3~3"D . In 
agreement with above consideration Ao (1Z = +1) is negative for any positive F . The negative 
Ai (1Z — — 1) appears when F is large enough. 

For sufficiently large F there will appear other solutions with negative A. However we restrict 
our consideration to the lowest solutions because they are most important for the problem of 



bound states. For illustration we present the eigenfunctions $(</?) for F = 1.4 and a = 7 in 
Fig. 5. The solution with 1Z — — 1 has two nodes (<£> = 0,7r), and in this sense it is "p-wave". 
There are no nodes in solution with TZ — +1, and in this sense it is "s-wave". 

Consider now the radial equation (|32|). For negative A there is the long-range attractive 
potential in Eq. (|32|) . As a result for every such A, A = Ao, Ai < 0, there appear the infinite 
series of bound states for the two holes. We can find their energies using the semiclassical 
approximation. 

E n « — l I leatp(-2--=(n + 3/4)). (37) 
Po v 4/IAI ' 



Here po is a short-range cut off parameter. The point is that the interaction — 1/p 2 is too 
singular at the origin and therefore it gives the "fall" of a particle to the center. To stabilize 
the solution we introduce infinite repulsion at p < p . It is obvious that the cut off parameter 
is due to the finite lattice spacing p ~ P1/P2 as well as to the contact interaction. To find 
the exact value of po one has to solve the equation (|i"7|) numerically with both short-range and 
long-range interactions taken into account. Below we will present examples of such calculations. 

Formula fl37|) is one of the most important results in the paper. It states that there are two 
(TZ = ±1) infinite series of bound states for two holes with S z — 0. However let us remember 
that there are two pockets in the Brillouin zone. Therefore we actually get four infinite sets 
of large-size bound states. The dependence of binding energies on the hole spin-wave coupling 
constant / and on the mass ratio a manyfests itself in Eq.([J7|) through the parameter A. For 
realistic value of t ~ 3 — 4 the constant F«l, see Fig. 2. It makes both Ao and Ai not large: 
I A 1 < 0.7, |Ai| < 1.5, see Fig. 4. Therefore the bound states (|37| ) are shallow and long-range 
that. 



Due to Eq.(|37j) typical momenta in these states are small: 

>(--^L(n + 3/4)). (3* 



1 1 
q ~ - ~ — exp l 

P Po V /\x\ 



Therefore the neglection of (E — e p — ejj in comparison with c<j q in Eq . fl2~T|) is justified. 



IV. SYMMETRY PROPERTIES 

Let us discuss the relation connecting the parity V and the quantum number TZ. The wave 
function of relative motion in coordinate representation is connected with that in momentum 
representation by relation fl2"5|). The inversed relation is 

2kOc5>(r n )e-* kr ". (39) 

n 

However let us remember that separation r n between two holes of opposite spins is always 
equal to the odd number of lattice spacings. The reason for this is that the ground state of the 
system is an antiferromagnetic one, see Eq.@. To avoid misunderstanding let us recall that we 
are considering dressed holes. Certainly due to the internal structure of the dressed hole fllOD 
there are configurations with the separation between two bare holes equal to the even number 
of lattice spacings. However r n is the separation between the centers of dressed holes, which is 
always odd. Thus only the terms with r n corresponding to the odd number of lattice spacings 
should be taken into account in sum (|3"9"D. Therefore 

9k- g = -#k, (40) 

where g = (±7T, ±7r) is the inverse vector of the magnetic lattice. Remember that the vector k 
is the deviation from ko which points to the band bottoms, ko = (±7r/2, ±7r/2). The inversion 
is the transformation k — > — k, ko — > — ko = ko — g. Therefore is transformed under the 
inversion as 

V : <7k - <7-k- g (41) 
Combining this with Eqs.([23|), ([!(]) we find 



V = -TZ (42) 

This is an important and surprising relation. Really, using the definition of 1Z (|23"D it is easy 
to verify that it may be presented as 1Z = ( — 1)', where I is the number of nodes of the angular 
wave function $(</?) in the region < ip < it. Eq.(|42|) states that the relation between I and V 
is unusual: V is positive (negative) for odd (even) 

It is useful to explane relation ( f4~2| ) using coordinate representation. The momentum with 
respect to the center of Brillouin zone is equal p = k + k, where k is deviation from the center 
of the pocket. Therefore 

* (r " ) = / 9kC ' pr W = e " k " r "/ 9ke ' kr "(&' (43) 

Due to Eq. (p5|) the last integral here is ip(r n ), and we get 

*(r n ) = exp(zk r n ) x ip(r n ) = exp(i^(n x + n y )^j x ip(r n ). (44) 



Here r n = (n x ,n y ). Since n x + n y is odd the exponent exp^|(n x . + n y )j is equal to ±z. It 
changes the sign under reflection of coordinate n x — > —n x , n y —n y . 

Relations d^0|) , (|42| ) permit one to find those representations of the symmetry of the system 
in which the bound states ( |3"TD should manifest themselves. The group of symmetry of the 
square lattice is C 4 „. It possesses four one- dimensional representations A±, A 2 , B ly B 2 of positive 



parity and one two-dimensional representation E of negative parity (see e.g. Ref. ||30|). 

Consider first the bound state with A = Ai. It has 1Z = — 1 and hence the positive parity 
V = 1. Therefore, we are to look for it among the one-dimensional representations. Further 
restriction comes from (f40p. In order to find it consider the nodes of the wave function g^. 
We know that the angular function $>\ 1 (ip) has two nodes, (p — 0, 7T, which are at the axis 
perpendicular to the face of the Brillouin zone, see Fig.5. The wave function in momentum 
representation has a similar angular behaviour. It vanishes when k is orthogonal to the face 



of the zone. This property does not contradict the symmetry in A 2 and Bi representations, see 
Fig. 6. For A\ and B 2 there is a trouble: the function satisfying (]40| ) in these representations 
must vanish at the edge, e.g. when k is along the boundary. This means that the series of 
bound states (|37D with A = Ai have to appear in A 2 and B\ representations only. 

A similar consideration shows that the bound states with A = Ao have to appear in two- 
dimensional E representation. 

We get the following full classification of the series of bound states (P7[): two of them with 
A = \x(TZ = —1,V = 1) are in B\ and A 2 representations, the other two with A = \q(TZ = 
1,V = — 1) are degenerate, they belong to E representation. 



V. NUMERICAL CALCULATIONS 

Now let us compare the analytical formula ([37]) with the exact numerical solutions of Eq.([T7|). 
Consider first the states of positive parity with TZ — — 1. One series of them belongs to the 
B\ representation. The symmetry of the wave function for this case is shown schematically in 
Fig. 6a. The short-range solution of this symmetry is usually called d-wave. The second solution 
belongs to the A 2 representation (Fig. 6b). 

The numerical solution of Eq. ( |l7|) is straightforward. We set the parameter a in the contact 
hole-hole interaction (|20"D , (|T2|) : a = 0.6; the mass ratio a = (3\/ (3 2 = 7. The value of hole-spin- 



wave coupling constant is taken from Fig. 2. The kernel T is calculated numerically with the 



help of Eqs.fllJD, (0), ©. 

A few first energy levels of B\ symmetry for different t are presented in Table 1. For each 
energy level we present also the corresponding value of (r 2 ) 

< r W^ d Witi 2 S- (45) 

For t = there is only one short-range bound state: = \^2(cosk x — cosk y ). It is due to 



the contact hole-hole interaction. In agreement with Refs. [24j-^7 the short-range bound state 
disappears at t ~ 2. 

The energies of large-size states are given in the second and the third sections of Table 1. 



The ratios E2/E3 are in agreement with formula E n /E n+ i = exp(2ir/ J\Xi\) which follows from 



Eq.(|37p. Numerical values of exp(27r/y |Ai|) with Ai determined by Figs. 4, 2 are presented in 
the last line of the Table 1. The state at t = 2 with E = —0.0417 is actually intermediate 
between short range and long range one. Therefore the agreement in ratio E 2 /E 3 is not so good 
as for the other states. 

There are no short range bound states of A 2 symmetry. The situation with long-range states 



is quite similar to that of B\ symmetry: they obey the Eq.(37). 

The only difference for the states of negative parity (TZ = +1) is that these states belong 
to the two-dimensional representation E, therefore each state is doubly degenerate. At t = 
there is only short-range bound state = 2 sin k x (or = 2 sin k y ), which is due to the contact 
hole-hole interaction. The short-range state (usually it is called p-wave) disappears at t « 2 



(see Refs. p4|-f27|j). At t 7^ there is the series of long-range states whose binding energies are 
described by Eq. flT7|) . 

We conclude that there is a good agreement between the results of direct numerical calcu- 
lations and the analytical approach developed in Sec. III. 

VI. CONCLUSION 

In the present work the long-range (I 3> lattice spacing) positive and negative parity two- 
hole bound states have been found for the holes on the quantum Neel background. The parity of 
these states is related to the angular wave function in a nonstandard way. Their small binding 
energy depends exponentially on the hole-spin-wave coupling. 



The natural question arises whether these states with the tiny binding energy could be 
relevant to high-T c superconductivity. There are the strong physical reasons which show that it 
is quite possible. We have demonstrated that there is the attraction between two holes caused 
by the spin- wave exchange. We have considered these two holes separately. If there is an 
ensemble of holes then one can expect that the attraction between a pair of two holes should be 
enhanced. The origin of this effect is the renormalization of spin-wave Green function caused 
by the particle-hole excitations. Therefore the binding energy of a pair is to be larger for the 
pair in the ensemble. The point is that a quite moderate enhancement of the interaction results 
in the very strong increase of the binding energy. 

To illustrate the later statement let us consider t = 3, and let us enhance the hole-spin-wave 
coupling constant by a factor of 1.5: / = 1.79 — > / = 2.69. Then solving Eq.(fL7|) with this 
value of coupling constant numerically we find the bound state of the pair belonging to B\ 
symmetry with the following characteristics: 

E = -0.103, (r?) = 8.4. (46) 

This is the ground state with a reasonable energy of E ~ 100f^ (we set J — 1 = 0.13eV). The 
state is stretched along one direction, and its maximal size is about 6-7 lattice spacings. 

The physical picture of the interaction enhancement is to be discussed in detail elsewhere. 
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TABLES 



TABLE I. 


The enerj 


;ies and the mean 


values of r 2 for 


a first few bound states of B\ 


symmetry 


t 





1 


1.5 


2 


3 


4 


E x 


-0.250 


-0.262 


-0.146 






~ 


(rl) 


1 


1.22 


1.95 








E2 


- 


-1.84- 10" 3 


-6.30 • 10~ 4 


-4.17- 10~ 2 


-2.21 • 10" 4 


-4.03 • 10" 6 


(rl) 




1.74 • 10 2 


6.34 • 10 2 


7.83 


2.65 • 10 3 


2 • 10 5 


E3 




—3.31 • 10 


—5.25 • 10~ 6 


-1.37- 10" 4 


w -2. • 10~ 7 


? 


(r 2 s) 




9.9 • 10 3 


8.5 • 10 4 


3.5 • 10 3 


w 2.8 • 10 6 


? 


(E2/ E^) num 




56 


120 


304 


1105 


? 






54 


123 


242 


1071 


2.15 • 10 4 



a The first section presents the small-size bound state which is due to the contact hole-hole 
interaction. It disappears at t ~ 2. 

b The large-size states are presented in the second and third sections. 

c In the last section the ratio (E 2 / E 3 ) numerical is compared with (E n / E n+l ) analytica i predicted by 
formula ([37]). 

d The accuracy of the code was insufficient to find extremely shallow levels denoted by "?" 



FIGURE CAPTIONS 

FIG. 1. The Brillouin zone of a hole in the t — J model. 



FIG. 2. The coupling constants. Solid line represents the hole-spin-wave coupling constant 
/ calculated in Ref. [29]. Dashed line represents the effective constant F = 2f 2 /(7i(3i) which 



governs the long-range interaction between two holes, see Eq.(p9|). 

FIG. 3. Single spin-wave exchange between two holes. Its contribution to the kernel of 
Eq . (|i"7|) is given in (|i~9|). 



Fig. 4. The eigenvalues A of the angular equations (0), ( |3~3"D versus F. The number at 
every curve is the value of the mass ratio a = /3i//?2=3,5,7,9. Fig. a - the solution Ao with 
1Z = +1. Fig.b - the solution Ai with 1Z = —1. Note that the vertical scale in Fig.b is different 
from that in Fig. a. 

Fig. 5. The angular wave functions for F = 2f 2 /(n(3i) = 1.4 and a = P1/P2 = 7. 

Solid line is the solution whose eigenvalue is Ai with 1Z = —1, and dashed line is the solution 
whose eigenvalue is Aq with TZ—1. 



Fig. 6. The symmetry of the positive parity (TZ = 
Bt type, Fig. b - A 2 type. 



— 1) bound state wave function. Fig. a - 



